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General  Remarks  about  the  Final  Report 

This  is  an  updated  and  enlarged  version  on  my  ONR-Report-6  sent  to  Professor  Pellat-Finet 
and  the  Program  Coordinator  on  October  1 1,  2002,  That  6*^  report  presented,  for  the  first  time, 
my  amplitude  matrix  approach  translated  into  quaternionic  form  My  matrix  approach  has 
been  concisely  described  in  the  earlier  elaborated  the  ONR-Report-5  (60  pages  text),  dated  31 
July  2002.  Its  comparatively  short  form  had  to  be  more  suitable  for  publication  in  the  capacity 
allowed.  Full  text  of  240  pages  has  been  presented  in  my  ONR-Report-3  (Final  Version)  dated 
1  June  2001.  That  was  the  Final  Report  of  the  previous  grant.  NOOO 14-00- 1-0620  (PR  # 
00PR06427-00  and  01PR07224-00), 

It  is  believed  that  the  text  of  this  Final  Report  fully  adopts  the  Kennaugh’s  inversion 
point  techniques  to  the  general  coherent  bi static  scatter  case.  Different  representations  of 
‘spinorial  scatterers’,  of  the  Jones  and  Sinclair  type,  on  the  Poincare  sphere  equipped  with  the 
Kennaugh’s  inversion  point  could  be  demonstrated  owing  to  use  of  the  time  and  spatial  frame 
reversal  concepts.  The  use  of  quaternions  (proposed  by  Pellat-Finet)  was  possible  taking 
advantage  of  their  close  relation  to  Pauli  matrices  by  which  amplitude  matrices  of  scatterers 
can  be  expressed.  To  stress  that  link  between  amplitude  matrices  and  their  quaternionic  forms 
the  same  symbols  have  been  used  for  them  with  the  only  difference  in  the  type  of  fonts:  the 
‘Times  New  Roman’  or  ‘Lucida  Handwritting’.  Therefore,  practically,  direct  translation  of  the 
matrix  to  quaternionic  formulae  was  possible  leading  to  the  end  result  in  the  form  strictly 
corresponding  to  those  with  the  Mueller  or  Kennaugh  matrices. 


Zbigniew  H.  Czyz 
Warsaw,  December  7,  2002 
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FOREWORD 

Fundamental  question  of  the  polarimetry  theory  in  communications  and  radar  is  how  to 
describe  polarization  of  antennas  and  waves,  the  process  of  polarization  transformation  when 
scattering,  and  how  to  express  the  received  voltage  with  its  phase  in  terms  of  polarization 
parameters  of  transmitting  and  receiving  antennas  and  those  of  scattering  object. 

Only  the  first  two  questions  have  been  answered  by  researchers  in  optics.  It  has  been 
established  that  the  polarization  and  phase  of  the  time-harmonic  locally  plane  waves 
propagating  in  a  vacuum  depend  on  three  real  parameters.  Two  of  them  describing 
polarization  can  be  presented  by  a  point  on  the  Poincare  sphere,  while  eight  real  parameters  of 
the  scattering  object  sufficiently  determine  transformation  of  polarization  and  phase  when 
scattering  ( [1],  [5]). 

The  problem  appeared  more  complex  in  the  case  of  partial  polarization  of  waves 
treated  as  the  incoherent  sum  of  those  completely  polarized.  Introduction  of  the  fourth  real 
parameter  appeared  necessary  together  with  16  real  parameters  of  the  scattering  object.  As  an 
example  of  such  presentation  can  serve  Stokes  4-vectors  and  Mueller  matrices  being  used  in  a 
statistical  optics  approach  ([1],  [2]). 

Next  problem  arose  in  communications  and  radar  considering  waves  propagating  in 
opposite  directions.  The  question  was,  for  example,  how  to  identify  three  real  parameters  of 
each  antenna  employed  in  the  direct  transmission  independently  of  its  destination  as  transmit 
or  receive  one,  and  how  to  use  those  six  parameters  to  express  amplitude  and  phase  of  the 
received  voltage.  One  of  satisfactory  solutions  to  this  problem  has  been  presented  in  this  text 
by  introduction  of  operation  of  the  local  spatial  frame  reversal  by  its  rotation  by  180®  about 
an  axis  perpendicular  to  the  propagation  axis.  Of  course,  such  approach  was  possible  when 
relating  polarization  to  the  spatial  frame.  That  can  be  considered  as  a  disadvantage  of  the 
method  but  the  results  obtained,  evidence  of  description  of  the  problem,  seem  to  justify  such 
an  approach. 

Having  main  idea,  mathematical  tools  have  to  be  chosen  leading  to  the  most  simple 
presentation  of  polarimetric  operations.  Those  tools  can  be  divided  into  two  categories; 
analytical  and  geometrical.  Geometrically,  the  elliptically  polarized  wave  can  be  represented 
by  a  polarization  helix,  here  proposed  and  better  serving  to  define  the  complete  polarization 
than  the  polarization  ellipse.  Kennaugh  [6]  extended  the  Poincare  sphere  concept  to  present 
geometrical  models  of  scattering  operators  by  introduction  the  concept  of  the  inversion  point 
located  inside  the  sphere.  That  concept  has  been  further  extended  to  broader  class  of 
scattering  operators  (originally  applied  to  backscattering)  what  has  been  reported  also  in  this 
text. 

There  are  many  ways  of  analytical  treatment  of  optical  and  radar  polarimetry. 
Interesting  reviews  have  been  presented  by  many  authors  dealing  with  the  radar  polarimetry 
theory.  In  opinion  of  this  author,  most  convenient  for  direct  computations  is  the  polarization 
and  phase  column  vectors  and  scattering  matrices  approach.  However,  there  exist  an  approach 
which  closely  links  polarimetry  with  relativistic  kinematics  described  by  Pellat-Finet  [10]. 
That  is  the  quaternion  approach  which,  what  should  be  stressed,  is  closely  related  to  the 
matrix  approach  (as  well  as  to  the  spinorial  approach  [7])  and  will  be  also  carefully  presented 
in  this  text.  Both  matrix  and  quatemionic  transformations  have  been  interpreted  geometrically 
by  using  models  of  scattering  operators  in  the  form  of  the  Poincar^  sphere  additionally 
equipped  with  the  inversion  point.  Also,  polarization  phasors  tangent  to  the  sphere  at  the 
polarization  points  appeared  necessary  to  precisely  determine  the  bases  of  polarization  and 
phase  vectors  or  quaternions. 
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1.  AN  INTRODUCTION.  THE  SPHERE  OF  TANGENTIAL  POLARIZATION 

PHASORS 

The  beneath  conducted  considerations  lead  to  construction  of  the  Poincare  sphere 
models  of  the  Jones  (Section  2)  and  Sinclair  (Section  3)  matrices,  or  their  corresponding 
quaternions,  presented  in  any  orthogonal  elliptical  basis  determined  by  one  tangential 
polarization  (TP)  phasor  B.  Its  orthogonal  mate,  the  collinear  Bx  phasor  of  the  same  phase, 
will  be  uniquely  determined. 

The  choice  of  some  special  bases  can  be  convenient  because  it  may  lead  to  simpler 
forms  of  transformation  matrices/quatemions  and  symmetrical  locations  of  special 
polarization  points  on  the  Poincare  sphere,  properly  rotated  versus  the  original  linear 
polarization  basis.  Also,  decomposition  of  the  Jones  and  Sinclair  matrices  into  product  of 
matrices  representing  real  and  imaginary  rotation  of  the  illuminating  polarization  and  phase 
(PP)  vector,  corresponding  to  the  TP  phasor  and  valid  in  any  such  a  basis,  enables  one  to 
indicate  intrinsic  properties  of  those  polarimetric  operators.  Therefore,  the  concept  of  the 
polarization  sphere  of  the  TP  phasors  and  its  phasor  basis  should  be  well  understood. 

The  TP  phasor  sphere  differs  from  the  known  Poincare  sphere  (of  polarization  points) 
by  introduction  of  additional  information  about  the  phase.  With  respect  to  the  polarization 
and  phase  space  of  the  PP  vectors  (compare  [8]),  the  TP  phasor  sphere  can  be  considered 

as  its  ‘squared  space’  with  all  three  angles,  determining  the  complete  polarization  and  phase, 
doubled. 

The  Poincare  sphere  and  the  sphere  of  the  TP  phasors  are  both  related  to  the  spatial 
frame  which  usually  is  the  local  right-handed  xyz  coordinate  system  with  the  z  coordinate 
chosen  along  the  propagation  axis.  The  equator  of  linear  polarization  points  on  the  Poincare 
sphere  of  unit  radius  corresponds  to  one  half  of  the  unit  circle  in  the  xy  plane.  On  the  sphere 
of  the  TP  phasors  one  deals  with  the  ‘double  equator’,  of  the  angular  length  of  4;r, 
corresponding  to  the  full  circle  in  the  xy  plane  and  resembling  the  Moebius  tape,  the  one¬ 
sided  surface.  The  TP  phasor  oriented  and  shifted  along  that  tape  comes  to  the  same  location 
at  the  tape  after  the  complete  4;r  shift.  Its  phase  doesn’t  change  when  shifting  along  the  tape, 
the  polarization  changes  only.  After  In  shift  it  arrives  at  the  ‘opposite  polarization’  point  but 
of  the  same  phase.  It  arrives  at  the  orthogonal  polarization  point  after  the  first  n  shift 
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Each  small  circle  of  that  sphere  is  also  of  the  4,n  lengths!  However,  when  shifted 
along  such  a  circle,  not  only  the  polarization  but  also  the  phase  is  changing.  The  sum  of  the 
double  change  of  polarization  (equal  to  the  solid  angle  subtended  by  the  circle,  in  radians) 
plus  the  double  change  of  phase  (the  effective  rotation  of  the  phasor  on  its  path)  equals  An . 
When  the  small  circle  reduces  to  a  point,  the  phase  is  only  that  what  changes  and  the  An 
rotation  represents  the  In  change  of  phase. 

Location  and  orientation  of  any  phasor  P  of  the  elliptical  polarization  (tangent  to  the 
Poincare  sphere  at  the  polarization  point  P)  can  be  determined  in  relation  to  another  phasor, 
for  example  the  original  basis  phasor  X  of  linear  polarization  and  null  phase  (tangent  to  the 
sphere  at  the  point  X  and  oriented  along  the  equator  of  linear  polarizations).  The  X  phasor 
represents  the  x  direction  of  polarization.  It  becomes  the  P  phasor  after: 

•  rotation  about  the  OX  axis  (called  also  the  axis)  of  the  sphere  by  the  angle 

of  the  right-hand-sense  (denoting  the  double  phase  delay  for  the  angle  of  positive 
value),  then 

•  shifting  in  direction  of  its  arrow  to  the  point  P  by  the  angle  of  ly^^r  (the  double 
angular  distance  between  polarizations  represented  by  the  X  and  P  points),  and 

•  the  rotation  about  the  OP  axis  (it  is  another  double  phase  delay). 


Ranges  of  those  angles  will  be  chosen  as  follows; 

Qi<2y  <n,  -n<2d<n,  -2n<2e<2n  (11) 

with  26-0  for  2y  =  0  and  2y  =  n .  The  couple  of  2y  and  25  angles  will  uniquely 
determine  any  polarization  point  on  the  Poincare  sphere.  Two  identically  oriented  phasors 
tangent  at  the  same  point  but  of  2e  parameter  differing  by  2n  will  be  considered  as 
representing  the  same  polarization  and  the  opposite  phase.  The  alternative  interpretation;  of 
the  same  phase  but  of  ‘opposite  polarization’  will  be  excluded  from  considerations  by  the 
above  proposed  ranges  of  angular  parameters. 

Any  basis  phasor  B  will  be  determined  by  other  three  angles,  2yl. ,  25^. ,  2e^. .  The 
second  phasor  of  that  basis,  Bx ,  will  be  uniquely  determined  by  shifting  B  to  its  antipodal 
point  on  the  Poincare  sphere  in  direction  indicated  by  the  B  phasor’s  arrow.  The  change  of 
the  TP  basis  X  ->  B  results  also  in  the  corresponding  change  of  the  Stokes  parameter  (right- 


Z.H.  Czyz  ONR-Report-7(FimI),  December  7,  2002,  NOOO 14-02-1-0222,  02PR04940-00  4 


handed  rectangular)  coordinate  system  Q^Ux  Yv  ^  basis,  the  Qg  axis 

will  cross  the  Poincare  sphere  at  the  B  point  and  the  new  equator  =  0,  as  well  as  the 
axis,  will  be  indicated  by  the  arrow  of  the  B  phasor. 

The  unit  column  polarization  and  phase  (PP)  vector  in  the  A"  basis  will  be  of  the  form 


cos  y  exp[-y(^  +  £•)] 

P 

r 

a 

sin^exp[-i-y(^-  £■)] 

X 

A. 

(1.2) 


The  change  of  the  TP  phasor  basis  X  ^  B  for  the  column  PP  vector  can  be  expressed  by  the 
unitary  matrix  transformation  (see  Appendix  A,  formula  (A. 2)): 


a 


h  «■ 


JB 


C^.  *u’’.  = 

A  "  A 


a' 


(1.3) 


<•;  detr;=+l 


—  h  a 

The  corresponding  column  Stokes  four-vector  depends  on  the  double  angular  parameters  ly 
and  2S  ; 


P.(^  = 


1 

p 

aa  *  +hh  * 

p 

"r 

1 

cos  2^ 

1 

aa  *  -hh  * 

1 

q 

sin  2y  cos  26 

ah  *  +ha * 

u 

sin  2y  sin  26 

A' 

j{ah  *  -ha*) 

A* 

V 

(1.4) 


Similar  change  of  the  TP  phasor  basis  can  be  expressed  by  the  following  orthogonal  matrix 
transformation  equation  (see  Appendix  E,  formula  (E.2)) 


p; ^p;  =D^p;  =D^.p;;  d^d^'  =D^Df.  =c/,«g{i,i,i,i} 
with  the  real  4x4  change-of-basis  matrix 


(1.5) 


m- 


1  0  0 
0  cos2;'  -  sin  2;^  cos  2£’ 

0  sin2;'cos2c5  cos2;'cos2^cos2£' -  sin  2^sin  2£' 
0  sin2/sin2(y  cos  2;^  sin  2(5  cos  2£’ +  cos  2(y  sin  2£’ 


10  0  0 

0  p  qR  qi 


0 

sin  ly  sin  le 

-  cos  2y  cos  2(5  sin  2e  -  sin  2S  cos  2e 

-  cos  2y  sin  26  sin  2e  +  cos  26  cos  2e 

(1.6) 
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with 

lp|=|qRl=|qil=l 

p*qR  =  qR  •qi  =  qi  •p  =  o>  (17) 

pxqR=qn  qRxqi  =  p»  qixp=qR 

and  where  is  a  vector  tangent  to  the  sphere  at  the  point  B  and  its  orientation  is  exactly  that 
of  the  B  phasor,  or  that  of  the  -  B  phasor  which  differs  from  B  by  the  angle  le  =  l7t . 

What  should  be  observed  is  that  the  orthogonal  transformation  matrix  depends  on  all 
three  double  angular  parameters,  though  the  four-vector  itself  bears  no  information  about  the 
phase.  That  is  the  result  of  rotation  of  the  Poincare  sphere  with  the  TP  phasors  of  the  PP  basis 
by  three  Euler  angles  2£^,,2y^,2S^  or  by  one  angle  2ip^  about  a  determined  n^.~*^  axis  of 
any  orientation  (for  details  see  again  Appendix  A,  substituting  phasors  X  for  P  and  B  for  K). 
For  that  reason  to  precisely  describe  the  orthogonal  elliptical  PP  basis  it  is  insufficient  to 
indicate  the  orthogonal  polarizations  only,  without  phases  of  the  PP  basis  vectors  determined 
by  directions  of  basis  phasors. 

2.  OPTICAL  POLARIMETRY.  THE  PROPAGATION  TRANSFORMATION 

2.1  The  Jones  matrix  decomposition  into  product  of  matrices 

The  time-harmonic  electromagnetic  polarimetry  considers  plane  waves  propagating  in 
a  homogeneous  isotropic  lossless  medium  in  which  they  can  be  sufficiently  described  by 
electric  vectors  only  because  their  magnetic  mates  are  unambiguously  determined.  Those 
waves  meet  obstacles  on  their  path  changing  their  amplitudes  polarizations  and  phases  being 
subjects  to  determination. 

Polarization  is  always  described  versus  an  established  local  spatial  frame,  a  right- 
handed  xyz  coordinate  system  with  the  z  axis  chosen  along  the  propagation  path.  In  optical 
polarimetry  z  axes  of  those  local  spatial  frames  are  always  directed  to  the  obstacle  for  the 
incoming  wave  and  out  of  the  obstacle  for  the  outgoing  (emerging)  wave.  So,  all  waves 
under  consideration  propagate  in  the  +z  directions  only.  Electric  vectors  of  those  waves, 
incoming  and  outgoing,  have  two  components,  along  x  andy  axes.  Those  components  used  to 
be  considered  as  complex  quantities  because  otherwise  it  would  be  impossible  to  separate  the 
three  interesting  terms:  the  real  amplitude,  ,  polarization  and  phase  (PP)  unit  complex 

vector,  //^  ,  and  the  wave  factor,  exp[y(G>?  -  kz)] .  Therefore,  the  real  components  are  treated 
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as  real  parts  of  complex  numbers  That  leads  to  the  following  expressions  for  electric  column 
vector  waves; 

r 

{I,z)  =  Kq{EIuI  QX'plJicot  - kz)]}  (2.1) 

Only  the  PP  vectors  are  subject  to  transformation  by  obstacles.  An  expression  of  such 
transformation  takes  the  form  of  an  equation  with  the  amplitude  2x2  complex  Jones  matrix 
(see[l],  [2],  [5]) 

{2.1) 


Here  XJ\So  are  tangential  phasors:  of  the  original  linear  polarization  basis  and  of  the 
incoming  and  outgoing  waves,  accordingly. 

Generally,  the  (amplitude)  Jones  matrix  any  R  basis,  will  be  presented  with 

two  indexes:  the  lower,  denoting  the  basis,  and  the  upper,  indicating  its  dependence  on  the 
Sinclair  matrix  A^  through  the  output  spatial  frame  reversal-by-rotation  matrix  Cg  : 


AO  -  (^O  *  A 


(see  Section  3  for  the  detailed  explanation). 


The  following  decomposition  of  the  Jones  matrix 


(2.3) 


(2.4) 


can  be  checked,  after  rather  tedious  work,  by  direct  substitution  of  the  following  rotation  and 
Lorentz  matrices  (see  also  Appendix  A,  formulae  (A. 3)  -  (A.7),  justifying  the  last  form  of  the 
rotation  matrix); 


fyROr° 

B 


[{A°e-^^°)+C^(A° 


Ale  +A( 
Ale  -  a; 
cos<^°  -  jn(  sin  (f)° 
0^3°  -^^sin  f 


*e^^°  A°e^^°-A° 

3  4 

A;e’^°+Al 

{-nl-jn\)sxn<f>° 
cos^°  +  y>/°  sin 


* 

*  e^^” 


B 


B 


(2.5) 
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with 


and 


cos^° 

Re{(7l°  +7l°)exp(-y^°)} 

0  •  jt  0 

Sin  (j) 

1 

-Im{(7l°-^nexp(-y^°)} 

n°  sin  (f>° 

-Im{(7i;  +7l3)exp(-y^°)} 

n°  sin 

B 

_Re{«-^3‘’)exp(-yr)} 

*A 


O 

B 


^0 

2\mI  +M°  +|dety4|  AlAl  *+A°A°  * 
~'^[a°A°*  +A° A°  *  M°  +  M°  +  |det  A\ 


1-Q  -U  +  yV 
-U-yV  1+Q 


(2.6) 


(2.7) 


where  Q^,Ul^  and  can  be  considered  as  coordinates  of  the  inversion  point  I  inside  the 

Poincare  sphere  of  unit  radius  in  the  B  basis  and  are  expressed  by  elements  of  the  first  raw  of 
the  Mueller  (and  Kennaugh)  matrix  as  follows  (see  Appendix  E): 


with 


b  =-L('M°  -M°  +M°  -M°  V  M°  =  A°  A°  *  )t  =  1234 

h  4-  /h  =  y4°  y4°  *4-y4°  y4°  * 

'-'3  ^7*^5  -^IB-^iB  ^-^AB-^IB 


and 


Coordinates 


(2.8) 

=M°,  +M°,  +m:,  +M°,  +  2 1  det  I  (2.9) 

and  are  called  the  normalized  Stokes  parameters  of  that  point 


There  is  another  way  to  decompose  the  Jones  matrix  expressed  in  any  orthogonal 
elliptical  basis  B,  by  a  matrix  Cg  of  the  ‘complex  pure  retarder’  of  unit  determinant 
multiplied  by  a  coefficient  (complex  number)  as  defined,  for  example,  in  Pellat-Finet  [10], 
Optik,  1992,  p.l06: 


with 


^B  = 


(2  10) 
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=1  det  Al  \ 

=arg4°  =  jargdet^^;  -n<^°<n 


(2  11) 


However,  in  general  case,  the  matrix  cannot  be  presented  in  such  a  form  This  is  because 
determinants  of  some  matrices  may  equal  zero  and  then  the  whole  matrix  vanishes 
identically.  Fortunately,  there  is  a  way  to  omit  that  difficulty.  Splitting  the  retarder  matrix  into 
product  of  pure  and  hyperbolic  rotation  matrices  (the  last  can  be  called  also  the  amplitude 
boost  matrix). 


with 


r  =r 


ROT° 


c 


BST  . 


det(’ 


=  detrf'  =1 


(2.12) 


cosh  a  +  m^g  sinh  a  (/Wjb  “  )  ^inh  a 

(m^g  +  jm^g )  sinh  a  cosh  a  -  m^g  sinh  a 


(2.13) 


one  can  introduce,  instead  of  ,  a  new  amplitude  matrix,  Af‘^ ,  owing  to  which  the  Jones 
matrix  will  not  vanish  identically  with  a„ .  Defining: 


one  obtains 


with 


and 


~  V  ^*0  '■  B  ~ 


1 


2  cosh  a 


C 


BST 


1  +  m^g  tanh  a 

(mig  -jm,g)  tanh  a 

2 

(m^B  +  .y/W3B)tanhQr 

1  -  m^g  tanh  a 

det  Af^  =  ^ =  -^(1  -  tanh'  a)  =  a^ 
4  cosh  a  4 

\SpattAf’^  =  ■“(!  +  tanh^  a)  =  jSpanAg  =  a, 

ag  =  SpanAg°‘^  +  2  det  A g°'^ 

=  2(a,  +  ao)  =  SpanAg  +  2  \  det  Ag  \ 


(2.14) 


(2  15) 


(2.16) 


(2.17) 


Finally,  the  Jones  matrix  can  be  presented  in  the  originally  proposed  form 
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A:=e^^°C 


ROT^  A  LOR 
B 


=  e^^ 


^ROT^ 


I  LOR 


(2.18) 


In  such  a  form  the  Jones  matrix  for  a,,  =|det/lg  |=0  will  not  vanish  identically  because  its 
coefficient  has  been  exchanged  for  /  2  which  is  always  greater  than  zero 

Further  (see  Section  2.7),  rQ=y[^j2  and  fo„=-Jo^/2  =  \  will  be  interpreted 

geometrically  as  radii  of  the  Poincare  sphere  models,  regular  and  normalized,  of  the  Jones  or 
Lorentz  matrices. 


2.2  Quaternionic  forms  of  the  Lorentz  and  rotation-after-Lorentz  amplitude 
transformation  matrices 


Quaternionic  expression  for  the  Jones  matrix  can  be  determined  by  introduction  of  the 
following  mutual  dependences  between  Pauli  matrices  and  their  quaternionic  versors  (see 
Pellat-Finet  [10],  Optik,  1992,  p.  101,  or  Misner  et  al  [7],  Chapter  41,  p.ll36): 


^0  = 


1  0 
0  1 


o  ^0  =  1,  a,  = 


1  0 
0  -1 


^1  = 


0  1 
1  0 


0  -j 
j  0 


Those  versors  undergo  multiplication  rules: 


(2  19) 


el^-\,  1,2,3;  =  ^3  =  ^  =  ^2  =  (2  20) 


In  terms  of  those  versors  the  quaternionic  form  of  the  Lorentz  amplitude  matrix  can  be  found 
by  applying  the  following  procedure  employing  formerly  established  dependence  between  the 
inversion  point  I  coordinates  and  components  of  the  mg  axis  of  the  imaginary  rotation 
(compare  the  last  formula  of  (2.7),  and  (2. 14)  with  (2. 17)): 


A 


LOR 

nB 


1-Q  -U  +  yV 
__U-yv  i  +  Q 

-(Q2:i  +Uo;2  +Va3  )J^ 
=a;o+(/n,o;,+/W2  2;2+'”3  2!3)B  tanha 


(2.21) 
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cosh  a 
1 

cosh  a 


cosh  a 


{cosh  «CTo+(»7,ct,+/W2  0^2 +"'3^3)5  sinh  a] 

(cos  /aCTo-/(/w,CT, +7^2 Si +'”3^3)5  sin  ja) 

exp(-ywgya)  =  — ^ — exp(wBa)  o  — ^ — exp{mgja) 


(2,22) 


cosh  a 


cosh  a 


with  the  matrix 


IEb  (QZi+Uct^+Vo;^))^  =  (qo;, +uo;2 +vo;2 ) g  jm^ 

-  tanh  a 


(2.23) 


The  /y  quaternion  corresponds  to  its  matrix  form  as  follows: 


=i-JiQ‘!,+Ve,+We,y^ 

=  1  +  J{m^e^  +  ^2^2  +  )b  a 
1 

exp(/WBya) 


(2.24) 


cosh  a 


with  the  (pure)  quaternion  m^  representing  the  axis  of  a  hyperbolic  (imaginary)  rotation  in 
the  B  basis 


mg  =  (/n,e,  +m^e^  +m^e^)g  s  (qe,  +0^2  +ve,)2 


(2.25) 


with 


q^  +u^  +  =  1 


(2.26) 


Those  normalized  Stokes  parameters,  coordinates  of  the  polarization  points  on  the  Poincare 
sphere  of  unit  radius,  can  be  expressed  in  terms  of  coordinates  of  the  M  point  representing 
polarization  of  the  incident  wave  corresponding  to  the  maximum  emerging  power; 


q 

M 

m, 

_  1 

■Q' 

u 

= 

m. 

—  1 

u 

L 

tanh  a 

V 

V 

B 

J^h_ 

B 

(2.27) 


Similarly,  the  amplitude  form  of  the  rotation  after  Lorentz  transformation  matrix 
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C 


ROT” 

B 


COS (f>°  -  jn°  sin  (f>°  (-/z®  -  yw®  )sin  (fp 

(n°  -  jnl)sm  (p°  cos ^®  +  jn°  sin 


=  cos (^° a Q- j{n^ a j+n° a 2'^n° a ^)g  sin 
=  cxp(-jnlf)  o  exp(n°J°) 


with  the  matrix  form  of  the  real  rotation  axis 


III  =  o;,  +n° +n° a,  )b  =  (qa,  +ua2  +va^  <=>  jnl 

can  serve  to  obtain  the  quatemionic  expression 

_  QQ5^o  +(n°e^  +«2^2  +'^3^3)b  sin 
=  exp(«°^®) 


(2.28) 


(2.29) 


(2.30) 


with  the  quaternion  n°g  of  the  real  rotation  axis,  through  an  A®  point  on  the  Poincare  sphere: 

"I  =  +«2^2  +"3^3)b  =  (q^l  +0^2  +  V£?3)f  (2.31) 


Of  course,  the  angle  of  rotation,  2(j>° ,  is  independent  of  the  PP  basis,  similarly  as 
tanh  a  =  (OI)„  does,  the  distance  from  the  inversion  point  I  to  the  center  O  of  the  Poincare 
sphere  of  unit  radius,  the  model  of  the  Jones  matrix. 


2.3  Quatemionic  form  of  the  Stokes  4-vector 

In  any  B  basis,  the  complex  quatemionic  form  of  the  unit  Stokes  4-vector  of  the 
completely  polarized  incident  wave,  represented  by  a  T  point  on  the  Poincare  sphere,  can  be 
found  through  the  density  matrix  expressed  by  the  PP  complex  vector  (see  [9]): 


=;^0  +  y(q^i  +ue^+ve,)l} 

'l  +  q  u-Jv' 
u  +  yv  1-q 

=  -^{0:0  +(qo;,  +ug:2  +vct3)I} 


JB 


(2.32) 


where  the  PP  column  vector  takes  the  form  (same  as  in  the  X  basis  -  see  Section  1). 
(Quaternions  proposed  earlier  by  Czyz  [3]  (1967)  did  not  use  their  complex  representation.) 
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cos^  exp[-y(^  +  s) 
sin  Y  exp[+ j{d  -  e) 


(2.33) 


A  digression.  The  unit  quaternion  of  the  Stokes  4-vector  for  the  partially  polarized  wave  will 
be  of  similar  form 

Pg  =  cos  2^  +  y(q  e^-^xxe^-^v  e^)\  sin  26  (2.34) 

corresponding  to  the  column  unit  Stokes  4-vector 


P  = 

*  n 


cos  2^ 

sin  26 coslYg 
sin  20  sin  2yg  cos2^j5 
sin  26  sin  2yg  sin  25 g 


-T 

cos  20 

T 

■  I  ■ 

sin  26  qg 

Q. 

sin  20  Ug 

u. 

sin  26  \  g 

1 

> 

_ J 

(2.35) 


with  tan  26  denoting  the  degree  of  polarization  which,  for  the  completely  polarized  wave,  is 
equal  to  tan  45®  =  1,  with  cos20  =  sin20  =-^.  Such  a  unit  4-vector,  or  the  corresponding 

quaternion,  can  be  represented  by  a  point  on  the  surface  of  a  hypersphere  of  unit  radius  in  the 
4-dimensional  space  of  Stokes  parameters.  The  crossection  of  that  hypersphere  by  a 
hyperplane  forms  the  Poincare  sphere  of  complete  polarizations  rotated  in  the  3- 

dimensional  space  according  to  the  implied  basis  B.  Such  representation  of  partial 
polarizations  by  points  in  the  4-dimensional  space,  instead  of  the  usual  representation  by 
points  inside  the  Poincare  sphere,  appears  useful  when  considering  reception  of  partially 
polarized  waves  expressed,  for  example,  by  dot  products  (see  Section'  )  of  the  corresponding 
quaternions. 

2.4  Quaternionic  forms  of  amplitude  rotation  matrices  changing  the  PP  vector  (the  TP 
phasor)  bases  of  vectors  and  matrices,  and  rotating  the  PP  vectors 

Quaternionic  form  of  equation  for  the  change  of  the  TP  phasor  basis  X  B  the  Stokes  4- 
vector  can  be  found  by  considering  at  first  the  matrix  in  the  X  basis  changing  the 
corresponding  PP  vector  (see  Appendix  A): 
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=[«f- 


X 

X 


(2.36) 


Quatemionic  forms  of  such  a  matrix,  rotating  the  PP  vector,  and  of  its  inverse  corresponding 
to  the  appropriate  matrices,  can  be  written  as  follows  (see  Appendices  A  and  D): 


Cx  =  cos<f>x  +(«,e,  +n^e^  sin<zJ® 

<=>Cx  =  cos (f>x  ao  -  j(P\  2:1  +  «2 5:2  +  "3 2:3  )x~"^  sin  <l>l 
Cg  —  Cx  *  =  cos^j-  —  (WjCj  +  /I2C2  ■*“  ^^3^3  );r  sin  (f)x 
<^Cg  =Cx*  =  cos^x^Q  +  +'^2[2  +”32;3)x^*  s\n<j)l 


(2.37) 


The  dependences 


cos^®  =  cos;k®  cos(<J  +  ff)® 
sin  <j)x  =  cos;k®  sin(<J  + 
sin  <f>x  =  -  sin  Yx  sin(<J  -  f)® 
sin  <j)l  =  sin  yl  cos(S  -  e)x 


(2.38) 


can  be  found  when  considering  the  equalities  (see  also  AppendixA): 


c‘^=h  »?]= 


a  -b* 

B 

cosy®  exp[-(<J  +  ]  -  sin  y®  exp[-(<J  -  e)x  ] 

_b  a* 

X 

sin  y ®  exp(<J  -  ^ )®  cosy^  exp(<J  +  ^ )® 

_ 


cos(f>x  -  sin  ^ 

/  ^^X >B  im.‘%X — \  — I*.  jL 


"x 

B 


(-/; 


X^B 

ZX 


■P^ix^ 

X^B 


)s\n(f>^ 


°  f  v«-«A  "-^B  iB 

cos  <Px  +  Jt^\x  sin  <Px 


(2.39) 

For  reasons  of  convenience,  the  rule  has  been  introduced  about  symbols  of  quaternions 
representing  transformation  matrices.  Symbols  of  mutually  corresponding  matrices  and 
quaternions  differ  only  by  one  defining  font.  For  example,  in  the  equivalence  of 

Cx^*  oCx  * ,  symbols  of  transposition  and  complex  conjugation  being  applied  to  matrices, 
together  with  indexes,  convey  without  change  to  quaternions  and  only  the  ‘Times  New 
Roman’  font  C  of  the  matrix  has  been  changed  for  the  ‘Lucida  Handwritting’  C  of  the 
quaternion  (both  in  italics). 


Matrix  equation  for  the  change  of  the  TP  phasor  basis 


V*.  jP  _ ^  X  ^.P  _ B  ik  P 

Ux  — ^  ~^X  X 


(2.40) 
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leads  to  the  following  equation  for  the  change  of  the  TP  phasor  basis  of  the  Stokes 
polarization  quaternion 

p;  -^p:=c,^  *p;c,^  cz>cl.  *ul.{r4yci.  (2.4i) 


Similarly,  the  equation  for  the  change  of  the  Stokes  polarization  quaternion  itself, 
corresponding  to  the  matrix  equation. 


^ClCl  *nl 


(2.42) 


IS 


fjP  v  SJV.  _  /^K  syV  ^  P  if  fpP^P  pK  if 

'B  '  B  ~  ^B  ^P  '  B  ^B  ~  ^B  ^B  'B'^B  '-'B 


(2.43) 


Quatemionic  equation  for  the  change  of  the  TP  phasor  basis  for  the  Jones  matrix,  also 
resembling  the  appropriate  matrix  equation,  is 


4° 


(2.44) 


2.5  Quatemionic  propagation  equation  in  the  form  of  the  product  of  quaternions  of  the 
Lorentz  and  rotation  transformations 


The  Jones  propagation  equation 


i'uf  =  Alul  ^Cr’ATul 


2  B 


(2.45) 


can  be  presented  in  the  quatemionic  form  through  the  following  mutual  relation: 


4 


where 


=  1  -  ;(Q^.  +  Ue,  +Ney 


(2.46) 


and 


(2.47) 
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0^9^  =  ^  {1  +  y(q  e,+ue^+\e,  } 

=  ^[1 


(249) 


exactly  corresponding  to  the  propagation  equation  cr^P*  =KgPj  with  the  Mueller  matrix 
K°  =  \J*{A°g  0i4^*)U  (see  Appendix  E). 


2.6  Derivation  of  the  Jones  matrix  decomposition,  in  the  quaternion  form,  using  the 
dephased  matrix  method 

The  dephased  Jones  matrix  in  the  quaternion  form  can  be  decomposed  into  product  of  two 
quaternions,  of  the  rotation-after-Lorentz  transformation  and  (not  normalized,  for  example) 


the  (pure)  Lorentz  transformation  (Lorentz  amplitude  transformation  without  real  rotation): 


(2.50) 


what  can  be  checked  by  expression  of  all  quatemionic  factors  through  elements  of  the  Jones 
matrix  and  by  accomplishment  of  quaternion  multiplications. 

To  perform  all  those  manipulations,  the  quatemionic  forms  will  be  useful  of  the  whole 
Jones  matrix  itself  and  its  versions:  transposed/conjugate,  transposed,  and  conjugate, 
accordingly,  which  are: 


(2.52) 
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+j{Al-A:)e,  +j{A;+A:)e,  +{Al  -  A:)e,}  (2.53) 

=  *+A°  *+J{A°  *-A°*)e,  +j{A°*+A:*)e,  -  iA°  * -A:*)e,}  (2.54) 

The  quaternionic  form  of  the  rotation-after-Lorentz  transformation  in  terms  of  the 
dephased  Jones  matrix  can  be  expressed  as  follows  (see  Appendix  B): 


/-Kor° 

'-'B 


{A°e-J^°  -e,A°  *e^  e^^° } 


(2.55) 


That  is  indeed  the  rotation-after-Lorentz  transformation  matrix  because  with 


ROT^  _ 


one  obtains 


y^ROT°  ROT°  Kc  _  /^  ROT°  y^ROT°  _ 


and  because  it  leads  to  the  quaternionic  form  of  the  Lorentz  transformation  matrix: 

^LOR  _  F'KOTo  if  -j^” 

r-lB  —  ^B  ^ 


(aj  +  a^  +y(bj^j 


2(a, +ao) 

f  2b,  2b3  2b,  1 

+  J 

'  c,  +  '  C2  j  -  e, 

o-o 

^  ^0  ^0  ^0  J 

B 

=  +u>,+v>3)J 


(2.56) 


(2.57) 


(258) 


and  to  its  normalized  version. 


A^g  —  {a|  +ao  +y(b]C,  +b3C2  yb5^3)i 

o-o 

=i-Xq;<’, +  u!,c’, +v>,). 


(2.59) 
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2.7  Poincare  sphere  geometrical  model  of  the  Jones  propagation  quaternion  (matrix) 

Simple  geometrical  constructions  present  the  Poincare  sphere  model  of  the  Jones 
propagation  matrix  (or  quaternion).  Two  special  functions  of  the  inversion  point  in  such  a 
model  should  be  observed: 

•  determination  of  the  emerging  intensity, 

•  explanation  of  the  polarization  and  phase  transformation. 

When  taking  the  diameter  of  the  sphere  equal  to  ,  the  outgoing  power  for  any 

polarization  point  T  of  the  incoming  wave  of  unit  power  equals  the  square  of  the  distance  of 
that  point  from  the  inversion  point  I: 

=(IT)^  =-^(IT)^  (2.60) 

4 

Such  a  result  can  be  obtained,  for  example,  when  considering  transformation  of  any  PP  unit 
vector  by  the  normalized  amplitude  Lorentz  matrix  in  an  M  basis  corresponding  to  the 
polarization  point  M  and  with  the  M  phasor  oriented  along  equator  in  the  plane  through  the 
points  I  and  T  in  direction  to  the  T  point: 


That  leads  to 


'l-Q 

0 

I 

a 

T 

■i  +  (oi)„ 

0 

cosy 

0 

1  +  Q 

u 

b 

M 

0 

l-(OI)„. 

sin  y 

[l+(OI)Jcos>-i 

,[l-(OI).]sin;-;: 


cosy 

n 

sin  y 

x;  =  Vi +(oi);  +2(01).  cos2r[,  =  ft 


or 

crj  =  1  +  (01)^  +  2(OI)„  COS  =  (ITy„ 


(2.61) 


(2.62) 

(2.63) 


what  agrees  with  the  Carnot  Theorem  for  the  triangle  OIT  inside  the  Poincare  sphere  of  unit 
radius  (for  (1M)„  =  1  +  (OI)„  and  2/^  being  an  angle  between  the  radii  (OT)„  and  (OM)„). 

Additionally,  the  result  of  the  above  normalized  Lorentz  transformation  determines  the 
following  relation  between  polarization  ratios  of  the  input  and  output  waves 


T.  1  -  (OI)„  T 


l+(OI)„ 


(2.64) 
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That  result  follows  also  similarity  of  triangles  NTI  and  MTT,  as  well  as  NIT’  and  MIT,  with 
the  TT’  chord  through  the  I  point.  It  resembles  the  Kennaugh’s  proof  of  similar  geometrical 
problem  for  the  inversion  matrix  (compare  his  M.Sc.  work,  p.  21,  eqn.  (75)).  In  this  case, 
however,  it  is  dealing  with  the  Lorentz  matrix.  It  allows  to  obtain  the  T"  =  T'x  point,  for  the  T 
point  chosen,  by  inversion  T  — >  T'  (through  the  point  I  back  to  the  sphere)  followed  by 
orthogonality  T'— >T",  geometrically  denoting  inversion  through  the  center  point  O  of  the 
sphere.  It  is  worth  noticing  that  the  Lorentz  transformation  does  not  change  the  phase. 

Similarity  of  triangles  NIT’  and  MIT  leads  also  to 


(IT')„  = 


l-(OI)^ 


(2.65) 


and  for  <f>"=  -  yj^ ,  with  the  T'T"  representing  diameter  of  the  sphere,  presents 

(IT')„+(IT)„  _  l-(OI)XAl)^ 


cos^"  = 


2Al 


(2.66) 


Altogether,  two  explanations  for  the  amplitude  Lorentz  transformation  exist:  the 
‘hyperbolic’  (or  equivalent  real)  T  ^T"  rotation,  or  the  successive  inversion  and 
orthogonality  transformation. 

According  to  its  second  function,  the  inversion  point  I  enables  one  to  find  the 
polarization  point  So  of  the  outgoing  wave.  To  obtain  that  result  one  can,  for  example, 

•  to  rotate  the  sphere  in  the  (previously  determined)  M  basis  about  the  axis  =  -1  by 


the  angle  2<f>"=2<f>l'  =  2ylj  -  2y]^ ,  dependent  on  the  T  phasor,  to  obtain  the  T"  phasor 
(after  the  Lorentz  transformation),  and  then 
•  to  rotate  the  sphere  about  the  tilf  axis  by  the  2<j>°  angle. 

It  is  interesting  to  add  that  also  the  last  axis  and  angle  of  rotation  can  be  obtained  having 
known  the  inversion  point  coordinates  in  a  special  basis,  characteristic  for  the  Jones  matrix. 
That  agrees  with  the  fact  that  the  dephased  Jones  matrix  depends  on  7  real  parameters  only. 
Four  of  them  are  the  Poincare  sphere  model  diameter  plus  three  coordinates  of  the  inversion 
point,  and  the  three  remaining  are,  for  example,  three  Euler  angles  by  which  the  sphere  should 
be  rotated  to  its  ‘characteristic’  orientation  (compare  similar  considerations  referred  to  the 
Sinclair  matrix  in  Czyz’s  ONR-Report-3,  Section  9.7  and  Appendices  G  -  J). 


To  arrive  at  the  So  phasor,  the  additional  rotation  of  such  rotated  phasor  is  required, 
by  the  -  2^°  angle  about  the  OSo  axis. 
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The  first  of  the  rotations  mentioned  above,  the  real  rotation  by  the  2(f>"  angle, 
corresponds  to  the  imaginary  rotation  determined  by  the  Lorentz  matrix.  It  depends  on  the  T 
point  chosen  and  replaces  the  amplitude  Lorentz  matrix  according  to  the  equalities: 


aLOR  _  . 

^nM  —  r,M  ’ 

=-l  '!<!>")  ■  (2.67) 

jm  T  _  T" 


The  last  transformation  can  be  rewritten  as  follows 


cos  ({>" 

sin  <f>” 

cos;" 

T 

cos  yl,  cos  ^^''+ sin  yl,  sin 

cos(/[,  -  f) 

cosyl' 

-  sin  (f>" 

cos,<f>" 

sin  y 

M 

sin  ylf  cos  </>"- cos  yl^  sin  (/>" 

sin(;^^  -f)_ 

_sin  rL\ 

(2.68) 

with  the  previously  obtained  angle  <f>" . 


3.  RADAR  POLARIMETRY.  THE  SCATTERING  TRANSFORMATION 

A  digression  about  the  E  (or  'amplitude  )  representation  of  the  electric  field  vector  (compare; 
Pellat-Finet  [8],  Optik,  1991,  pp.27,  and  70). 


The  electric  vector’s  plane  wave  propagating  in  the  +  z  direction  of  the  local  right- 
handed  xyz  spatial  frame  can  be  written  in  different  forms,  depending  on  the  choice  of  the 
orthogonal  basis  of  collinear  phasors,  X  or  B,  in  which  its  angular  parameters  are  being 
expressed; 

^E^{(z)  =  El[l^  exp[y(ai/-fc)]} 

=  ^0  [lx  1, ]Re{(cos;-^;/^  +  sin  y^uf  exp{j2S^ ))exp[j{cot  -{kz  +  v^ ))]} 

(3.1) 

or  in  terms  of  the  B  basis  column  PP  vectors 


with 


X  _ 
V  “ 


/^]Re{C>^  exp[7(arf  -  kz)]} 

/^]Re{(cosrB?/x  +  sin  yluf  exp{j2dl))exp[j{oyt  -  {kz  +  ))]} 

(3.2) 


T 

.  Xk 

'o' 

0_ 

,  u^  = 

1 

p 

X  > 


=  S^+s 


p 

B 


(33) 
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(3.4) 


(3.5) 


(3.6) 


The  same  wave,  when  changing  the  propagation  direction  for  the  opposite  one  {and  when 
applying  the  'time  reversal'  concept)  can  be  successively  written  as  follows: 


(t,z)  =  El[l^  i^]Re{7/^  exp[y(-<yt-/rz)]}  =  £o  /^]Re{(?/^.)*  exp[y(<yt +  ^z)]} 
=  ^0  [k  Jy  ]Re{(cos/f-?/;^  +  sin  yl-uf  exp(-y2(J^ ))exp[y(rtV  +  {kz  +  ))]} 


(3.7) 


or  in  terms  of  the  B  basis  conjugate  column  PP  vectors: 


==  El  [/,  ]Re{(cos  rl  {ul- )  *  +  sin  yl  (?/f )  *  exp(-y 2^^^ ))  exp[y(m/  +  {kz  +  ))] } 

(3.8) 


Both  waves  have  the  same  polarization,  represented  by  the  P  point  on  the  Poincare  sphere, 
and  their  equations  present  two  polarization  helices,  coinciding  at  the  moment  of  time  t  =  0 
but  moving  in  opposite  directions.  Position  of  such  a  helix  in  space  at  the  time  t  =  0, 
expressed  versus  the  z  axis  (the  propagation  axis),  can  be  called  its  spatial  phase ,  for  example 
or  Vg,  eventually  el  or  e^ ,  depending  on  the  definition  applied  and  on  the  basis 
assumed.  However,  senses  of  polarization  ellipses  traced  by  both  waves  penetrating  the  2  =  0 
plane  are  different!  Therefore  it  is  not  advisable,  in  radar  applicatiotts,  to  define  the 
polarization  by  the  polarization  ellipse.  Instead,  the  polarization  helix  together  with  its  spatial 
phase,  both  determined  by  the  same  TP  phasor  P  on  the  polarization  sphere,  is  the  best 
solution  to  that  problem.  In  time,  such  a  phasor  rotates  in  opposite  directions  depending  on  the 
direction  of  propagation,  ±  z 
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3.1  Spatial  frame  reversal-by-rotation  and  orthogonality  transformations 

The  Sinclair  (‘scattering’)  matrix  can  be  obtained  from  the  Jones  (‘propagation’)  matrix 
through  the  reversal-by-rotation  of  the  output  local  spatial  frame.  That  results  in  the  output 
complex  amplitude  (CA)  unit  vector  represented  by  the  conjugate  PP  vector,  wf  *,  of  the 
scattered  wave; 

Ay,=X^ul*-  A,=C°A°  (3.9) 

In  the  case  of  the  Sinclair  matrix,  the  emerging  wave  propagates  in  the  -  z  direction  of  the 
local  spatial  frame,  ‘reversed -by-rotation’  versus  that  applied  for  the  Jones  matrix. 

The  transformation  equation  for  the  spatial  frame  reversal  by  the  180®  rotation  about 
the  axis  is  of  similar  form  as  that  with  the  Ag  matrix  (see  also  Appendix  A): 

w^  =(a^ =  (ab  * +ba*)^ 

For  the  matrix,  its  parameters  are;  w=l,  and  u=0. 

Other  special  properties  of  such  a  ‘Sinclair  matrix’,  Cg,  which  can  be  called  the  ‘free 
space  scattering  matrix’,  are;  its  symmetry,  and  =1,  because  the  inversion  point  is  in  the 

center  of  its  normalized  Poincare  sphere  model.  That  is  the  reason  for  applying  the 
symbol  for  that  matrix. 

Both  matrices  transform  identically  under  the  change  of  the  PP  basis; 

Ag^iAyi,  cycyici  (3.11) 

what  can  be  explained  by  inspection  of  the  transmission  equation; 

Vr  =uUgul  ={uy^g){Cl.AgC-^){ClMl)  =  u^Ay,  (3.12) 


D 


detC°  =-l. 


(3.10) 


Only  the  orthogonality  transformation,  also  of  the  same  form. 
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C 


C’ 


0  -1 

1  0 


(3.13) 


does  not  depend  on  the  PP  basis  because 


>< 

II 

>< 

(3.14) 

The  following  quatemionic  expressions  can  be  found:  for  the  reversal-by-rotation  matrix,  for 

its  conjugate  form,  and  for  the  orthogonality  matrix. 

Cg  -  y(-  Im  w  -  Re  w  e,  +  ue^ 

(3.15) 

Cg  *  =  y(lm  w  -  Re  w  e,  +  ue^ 

(3.16) 

=e, 

(3.17) 

In  scattering,  the  transformation  equations  in  the  quatemionic  form  should  take  into 
account  conjugate  values  of  some  PP  vectors.  Therefore,  to  arrive  at  the  Stokes  polarization 
quaternion  in  the  reversed-by-rotation  spatial  frame,  one  should  write 

V2  »»«/  •  =  C|  •  -Jluf  o  =  c;  •  (3.18) 

where: 

=7j{l  +  y(qe,  +\ie^  +ve3)|}  (3.19) 

and 

=  TT  +  M  e,+ue,-y  e,  }  (3.20) 

Those  quaternions  correspond  to  the  complex  amplitude  (CA)  vectors  of  waves  propagating 
in  the  +z  and  -z  directions  of  the  local  spatial  frame,  accordingly.  Only  the  ‘CA  quaternion’ 
of  a  wave  propagating  in  the  +z  direction  equals  the  Stokes  polarization  quaternion, 

Pb"*  (3.21) 

Therefore,  in  the  case  of  the  ‘scattered’  emerging  waves  which  propagate  in  the  -z  direction 

of  the  local  spatial  frame  (contrary  to  the  ‘propagated’  emerging  waves  represented  by  the  CA 
and  PP  vectors)  their  Stokes  polarization  quaternions  correspond  to  the  CA  vectors,  //g. 
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for  the  hypothetically  reversed  direction  of  propagation  of  those  waves.  After  all  quatemionic 
multiplications  one  obtains  the  end  result. 


Pb' J^i  =;^{l  +  y(qe, +uej +ve3)|}  (3.22) 

where 


pSo 


T 

So 

"1 

0 

0 

0 

q 

and  D®  = 

0 

l-2u^ 

-  2u  Re  w 

u  ^  -  Re  w  ^ 

-  2u  Im  w 

u 

0 

-  2u  Re  w 

-  Im  w^ 

V 

B 

_0 

-  2u  Im  w 

-  Im  w^ 

u^  +  Rew^ 

(3.23) 


Similarly,  the  orthogonality  transformation  is 


pj-  =c^p/-c^ 


—  Cj  ^  {1  +  y(q  Cj  +  u  ^2  ^3) 


~  ^  y(Q  ^1  +  n  Cj  ■^'^^3)5) 


(3.24) 


3.2  Derivation  of  the  Sinclair  matrix  decomposition  in  the  quaternion  form  using  the 
dephased  matrix  method 

The  quatemionic  form  of  the  rotation-after-inversion  transformation  matrix  in  terms  of 
the  dephased  Sinclair  matrix  (compare  Appendix  B)  is: 


with  its  inverse 


y^ROT  «  _ 
*-'B  - 


^3  +^3^  } 


a 


ROT 


*eJ^ 


-4^3  ) 


(3.25) 


(3.26) 


This  is  indeed  the  rotation-after-inversion  transformation  because 
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and  because  it  leads  to  the  quatemionic  form  of  the  inversion  transformation  which 
corresponds  to  the  Lorentz  transformation  followed  by  the  inversion  through  the  center  of  the 
Poincare  sphere,  the  inverse  versus  the  orthogonality  transformation: 


=  -e. 


*-e^  Ag  Cj  e-^^^}Ag 


(3.28) 


When  searching  the  derivation  of  the  Lorentz  matrix  observe  please  that  the  form  of  the 
inversion  matrix,  or  its  quaternion,  must  not  depend  on  the  choice  of  the  Sinclair  or  Jones 
matrices  taken  together  with  their  phase  terms.  So: 


Ar  = 


-e. 


{Ag°  *  -e.  A, 


'B  ^  f'  'B 


-e,Ag 


LOR 


(3.29) 


This  is  because  (see  also  Section  3.4): 


{Ag°  *-e,Ag°e,  e' }Ag°  ^  a,  +  +  j{h,e,  +  h,e^  -  Jb,e,)g  =  {Ag  *-e,  Ag  e, 


(3.30) 


The  normalized  version  of  the  inversion  matrix  quaternion  becomes 


A„g  —  Cj  (a,  +  a^  +y(b|ej  +b3e2 
<^0 

+u;e,+v>,).) 


(3.31) 


LOR 


3.3  The  dependence  of  the  inversion  and  rotation-after-inversion  matrices  on  elements 
of  the  Sinclair  matrix 

In  order  to  find  mutual  dependences  between  elements  of  the  Sinclair  matrix  and  the 
axis  and  angle  of  rotation  after  the  inversion  transformation,  it  is  instructive  to  decompose  that 
matrix  as  follows  (see  Czyz’s  ONR-Report-5,  July  31,  2002,  p.26,  Eqs.  (3.60-61)): 


Z.H.  Czyz  ONR-Report-7(Final),  December  7,  2002,  NOOO 1 4-02- 1-0222,  02PR04940-00  25 


with 


As- 


"^3 

L^4 


^qROT 


ilNV 


Cr  •  =  +  C'  (a,  •  e'< )} 


(3.32) 


0  L 


^  ^  ^ 


A^e 


-ji 


+  ^ 


* 


,-Ji 

-  A^e~^^  +  A-,*  e’^ 


JB 


COS  (f>  +  7’//,  sin  (f>  (-//j  +  y'//^ )  sin  (f> 
(//j  +  jn^ )  sin  (j)  cos  (f>  -  jn^  sin  (j) 


(3.33) 


and 


A'”’  =e  •+A,C^e-“<}4 

/CT, 


CTn 


cr. 


A2A^  *+^4/ii  *  A/3  +A/i  +  |det^ 

-M2 -M,-\i3itXA\  -A,A2  *-A,A 


* 

4  JB 


ba+yl),  ^  +  b, 


--^-b,  -b3+7l), 


-U-7V  1  +  Q 

-1  +  Q  U-7V 


3 

nl 


(3.34) 


-  nB 


By  comparing  the  two  last  forms  of  the  amplitude  matrix  of  rotation-after-inversion,  one 
obtains: 


(3.35) 


cos^ 

Re{(^3  -A^)b  exp(-7^)} 

//,g  sin  (j) 

1 

Im{(^3  +A^)b  exp(-7^)} 

n,B  sin  (f> 

Im{(^,  -4)Bexp(-7^)} 

Jhs  sin  (!>_ 

Re((^,  +^2)b  exp(-7^)}_ 

Of  course,  the  angle  of  rotation  after  inversion,  2^,  does  not  depend  on  the  B  basis.  This  is 
because  the  difference  of  elements  of  the  Sinclair  matrix  second  diagonal  is  independent  of 
that  PP  basis. 

To  summarize  mutual  dependences  between  matrices  of  decompositions  of  the  Jones 
and  Sinclair  matrices,  it  is  instructive  to  •write  (see  Appendix  B,  formulae  (B.3)  and  (B  5)): 
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r’f  ^ *  =  -jCH'i  * *  =  -jcicf^°c 


X  a^OR 


A'f  =-r’'/4g 


(3.36) 


That  is  an  immediate  result  of  the  equalities: 


(3,37) 


with 


e  -  p°  j.  i 


(3,38) 


3.4  Similarities  and  differences  of  the  Poincare  sphere  models  of  the  Sinclair  and  Jones 

matrices 

For 


and 


Al 


^2 

a;  A,° 


=  Q*^b 


As 


C°A 


o 

B 


(3.39) 

(3.40) 


^  =  1,2, 3,4 


(3.41) 


it  can  be  checked  immediately  that  (see  also  Appendix  E): 


Mjg  +  Mag  “  ^  2B  "^^43  “^1 
Mjg  +  M,g  — Mjg  +  Mjg  — 3]  “b,g 

A°2bA°s  * +A:sA(s*  =  A,,A,,  *+A,,A,,*  =  h,,  +Jh 
I  A°sA(s  -  A°,A:,  1=1  A,, A,,  -  A,,  A,,  |=|  det  A°  |=|  det  A 


B  I 


(3.42) 


what  results  in 


i(M°g  +M°g  +M:g  +M°g)  =  i(M,g  +M3g  +Mag  +M,g)  =  a, 
i(M°g  -  M°g  +  M:^  -  M^g)  -  i(M,g  -  M3g  +  Mag  -  M,g)  =  b,g 


(3.43) 
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and  means  that  the  Jones  and  Sinclair  matrices  have  the  Poincare  sphere  models  of  the  same 
diameter, 


=  72(a,  +ao) 


(3.44) 


and  the  same  coordinates  of  the  I  and  M  points  in  the  same  B  basis  (compare  (2.27)): 


r 

b, 

'q' 

-2 

u 

=  -  tanh  a 

u 

_v_ 

0 

nB 

_b5. 

B 

V 

(3,45) 


with  the  same  independent  of  the  B  basis  parameter 

tanha  =  2^h], ! a,  ^IhJ a,  =  (01)„  (3.46) 

Results  obtained  are  physically  obvious  because  emerging  power  cannot  depend  on 
the  spatial  frame  reversal  for  the  outgoing  wave.  Also  they  indicate  the  a  angle  to  be  a 
common,  intrinsic  parameter  of  the  Jones  and  Sinclair  matrices,  uniquely  determining  the 
inversion  I  point  coordinates  of  which  depend  on  the  PP  basis  B,  arbitrarily  chosen  That 
inversion  point  determines,  in  a  different  way,  the  inversion  and  the  Lorentz 
matrices/quaternions  being  used  to  build  the  Sinclair  and  Jones  matrices/quaternions, 
respectively. 

Contrary  to  that,  the  rotation-after  inversion  and  the  rotation-after-Lorentz 
transformation  matrices/quatemions  are  different.  Their  axes  and  angles  of  rotation  depend  on 
the  reversal-by-rotation  transformation  of  the  spatial  frame  for  the  outgoing  wave.  General 
dependences,  for  any  PP  basis  5,  joining  axes  and  angles  /?g,2^  with  n°g,l<f)°  are  as  follows; 

cos^° 

«,  sin  <(> 

«2  sin  <f)° 
n°  sin 


■  n^g  Re  w_j. 


Im  w^ 


u^  cos  <f)^{n^g  Rew^.  -n^g 


iB  ““  v,^)sin^ 

Im  w^  )sin  <f) 


Re  w^  cos^  +  (w,g  Im  w^  -  n^gU  ,^ 


Im  w^  cos<f) 


(«2bU^  -n,g  Rew^ 


^  )  sin  ^ 
^  )sin  ^ 


(3.47) 


They  directly  correspond  to  the  quatemionic  form  of  the  spatial  frame  reversal-by-rotation 


matrix 
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cr  -jC^QC, 


ROT 
B 


-  Im  w  -  Rew  +  ut?2)v  ][cos^  +  +  n.^e^  +  n^e^  )g  sin 

=  (ue,  +  Re  w  ^2  +  Im  w  e3)f.[cos^  +  (//,e,  +  n^e^  +  sin  (j)'\ 


(3.48) 


=  -(//,suf-  +  n-^B  Re  wf.  +  Im  wf-  )sin  (j) 

+  e,[uf,  cos^  +  (//3g  Rewf.  -//2s  Imwf.  )sin  (p] 
+  e2[Re  wf.  cos^  +  (//,a  Im  wf-  -  //33U^- )sin  (pi] 
+  e3[Im  wf.  cos^  +  -  n^g  Re  w )sin 


(3.49) 


=  cos^°  +{n°ge^  +n°ge^  +/;33e3,)sin 


In  the  special  case  B  X (in  the  linear  PP  basis  X,  where  w=l  and  u=0)  one  obtains. 


cos^° 

-  /;2  V  sin  (p 

-Im[(/1,  -A^)^.e 

Re[(/l2°+/ir),.c' 

n°  sin  (P° 

//3  V  sin  (p 

1 

Re[(^,  +  A,) j^.e 

1 

-Im[(.42°-4°).v‘-’ 

n\  sin  (p° 

COS(p 

Kc[{A,-A,),.e 

-Im[(/i:  +/I0.vt’ 

n°  sin  (p° 

B 

.s 

w, 

1 

_ 1 

Re[(.4:-.43°),c^^‘°]_ 

(3.50) 


3.5  Reception  of  the  scattered  power 

The  power  received  by  an  antenna  from  a  scattered  wave  can  be  expressed  by  the  ‘dot 
product’  of  two  Stokes  polarization  quaternions,  one  of  the  antenna  and  another  one  of  the 
scattered  wave.  Such  expression  follows  its  corresponding  matrix  form: 

=  P^Pl  =i(l  +  q^q^  +u^u^  +v^^)  (3.51) 

In  case  of  the  partially  polarized  scattered  wave  the  received  power  expressed  by  the  product 
of  the  Stokes  polarization  vectors/quaternions  is  (when  applying  forms  (1.4)  and  (2.35)): 


/^  =  P."Pb  =i^{cos20^  +{qWB+«  +  v^v^)sin  2^" }  (3.51a) 
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4.  DIRECT  TRANSMISSION  BETWEEN  ANTENNAS 

When  considering  the  direct  transmission  between  antennas  it  seems  natural  to  express  the  PP 
vectors  of  those  antennas,  or  their  corresponding  Stokes  polarization  quaternions,  in  their 
‘own’  local  spatial  frames  xyz,  with  the  z  axes  directed  outward,  ‘from  antennas’  looking  at 
each  other,  and  with  the_y  axes  remaining  parallel  (the  last  condition  is  suitable  for  applying 
the  C°j(  matrix,  or  the  quaternion,  in  the  here  assumed  form).  Following  such 
assumptions,  the  received  power  can  be  presented  in  the  form: 


Pr  =\KClnl  \'  =  Ps^lK  =PbPb“  =PrPl 
- 

~'B  'B  ^B 


(4.1) 


Here,  the  polarization  quaternion  of  one  antenna,  for  example  corresponding  to  the  Ug  as  its 
unit  complex  height  or  its  polarization  and  phase  vector,  is 


p." 


1  +  q  u-yv 
u  +  yv  1  -  q  _ 


-\R 


(4.2) 


The  Stokes  quaternion  corresponding  to  the  conjugate  complex  height  of  another  antenna. 


PJ-  <z>42rtl*ul  = 


1  +  q  u  +  yv 
u-yv  1-q 


=;^(l  +  y(q^i  +uej  -\e,Yg) 


JB 


(4.3) 

(which  is  not  a  polarization  quaternion),  determines  the  polarization  quaternion  of  that 
antenna  when  expressed  in  the  spatial  coordinate  system  reversed,  by  rotation,  versus  its 
‘own’  local  spatial  frame. 

uTo  „  ro~ro*  1 

Pg  <Z>V2»BWfl 

_  /^o  if  «r-^o 
~^B  '  B  ^B 

-  y(Im  w  -  Re  w  e,  +  uCj  )^  ^  (1  +  y(q  e,  +  u  -ve^Yg}  y(-  Im  w  -  Re  w  e,  +  uCj  )f. 


1  +  q  u-yv 
u+yv  1-q 


^  “ --^(l  +  y(qe, +UC2 +ve3)g°} 


-  B 


=  ^  (1  +  ye, [(1  -  2u ^  )^  q^  -  (2u  Re  w)> ^  -  (2u  Im  w)^-  v",  ] 

+  yej[-(2uRew)^q^  +(u^  -Rew^)®  u^  -(Im  w^xv^] 

+  ye3[-(2u  Im  w)®  q^  -  (Im  w^)®.  u^  +  (u^  +  Re  w^)®  v^] 

(4.4) 
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Observe,  please,  that  physically  the  same  polarization,  emitted  by  the  same  antenna  in  the 
same  direction,  can  be  presented  by  two  different  polarization  4-vectors  (quaternions)  when 
expressed  in  mutually  reversed  (by  rotation)  spatial  coordinate  systems.  It  confirms  the  rule 
that  description  of  polarization  of  an  antenna  or  wave  is  always  referred  to  its  spatial  frame. 

For  the  linear  X  basis  of  collinear  phasors,  or  for  the  left  circular  (LC,  RC)  basis  of 
parallel  phasors,  for  which  w  =  1  and  u  =  0,  one  obtains  the  known  result, 

P.  =(^'  =(.P“°  •P’)x,ivc.nc, 

(About  the  parallel  phasor  basis  see:  the  Czyz’s  paper:  “A  comparison  of  different  kinds  of 
orthogonal  polarimetric  bases”  in  MIKON  2002  Conf  Proc.,  pp.  164-167,  or  his  ONR-Report- 
3,  pp.53-54.) 

It  should  be  stressed  that  the  ‘dot  product’  of  polarization  quaternions  in  the  above 
presented  formulae  is  different  from  that  proposed  in  [9]  or  [10]  because  it  is  identical  to  the 
Euklidean  scalar  product  of  Stokes  four-vectors  or  to  the  real  part  (in  quatemionic  sense)  of 
the  complex  quaternions’  product.  So,  in  applications  to  transmission  between  antennas  or  to 
reception  of  waves,  complex  quaternions  corresponding  to  the  Stokes  4-vectors  are  not 
elements  of  the  Minkowski  space  with  Lorentz  metric.  They  undergo  the  Cartesian  metric, 
also  for  hypercone  of  complete  polarizations  resembling  the  light  cone  in  the  theory  of  special 
relativity.  Nevertheless,  they  undergo  Lorentz  transformations  discussed  in  [10]  as  vectors  in 
relativistic  kinematics. 

5.  GEOMETRICAL  MODELS  OF  THE  PROPAGATION  AND  SCATTERING 
AMPLITUDE  MATRICES  AND  THEIR  COMPARISON 

Models  of  matrices  are  tangential  polarization  phasor  (TPP)  spheres  of  given  diameter, 
inversion  point,  axis  of  rotation,  angle  of  rotation  about  that  axis,  first  phasor,  B  say,  of  the 
orthogonal  collinear  basis,  and  angle  of  the  ‘general’  phase.  Common  to  both  propagation 
(Jones,  Ag )  and  scattering  (Sinclair,  Ag )  matrices,  corresponding  to  the  same  obstacle  being 

met  by  the  wave  on  its  path,  is  diameter  of  the  sphere,  ,  the  inversion  point,  I,  and 

(usually)  the  basis  assumed  in  which  elements  of  matrices  and  three  coordinates  of  the  I  point 
are  determined.  Different  are:  the  way  of  using  the  inversion  point,  and  the  axes  and  angles  of 
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rotation.  General  phase  differs  by  90°  for  the  basis  assumed.  First  transformation  of  the 
incident  wave’s  TP  phasor  is  with  the  use  of  the  I  point,  then  comes  rotation. 

The  simplest  is  the  inversion  transformation,  being  applied  when  considering 
scattering,  determined  by  the  inversion  matrix,  A^g'' .  It  moves  the  polarization  point  T  of  the 
incident  wave  by  projection  through  the  inversion  point  I  back  to  the  sphere,  to  the  point  T’=T 
INV.  Direction  of  the  transformed  phasor,  tangent  at  that  point,  can  be  determined  by  the 
following  procedure.  The  two  points,  T  and  T’,  should  be  joined  by  the  oriented  shorter  great 
circle  arc  in  such  a  way  that  its  origin  is  placed  at  the  inverted  point,  T’,  and  its  end  at  the 
incident  polarization  point,  T.  At  first,  the  incident  phasor  should  be  rotated  to  become 
collinear  and  identically  oriented  with  the  T’T  arc,  then  it  should  be  shifted  ‘back’  to  the  point 
T’  and  rotated  again  in  the  same  direction  by  the  same  angle.  After  such  procedure  it  will  be 
oriented  at  the  same  angle  but  of  opposite  sign  versus  the  T’T  arc  in  comparison  with  its 
previous  direction.  Analytically,  such  result  follows  the  conjugate  value  of  the  inverted  PP 
vector  equal  to  CA  of  the  transformed  wave.  At  the  end  the  phasor  should  be  multiplied  by 
the  real  number  equal  to  the  IT  distance  inside  the  sphere  of  diameter  . 

More  complex  is  the  Lorentz  transformation  being  applied  when  considering  the 
propagation  transformation  and  determined  by  the  Lorentz  matrix,  .  At  first  it  operates 
like  the  inversion  matrix  but  afterward  it  realizes  the  orthogonality  which  is  the  reversed 
transformation  versus  inversion  through  the  center  point  0  of  the  sphere.  If  the  I  and  0  points 
coincide  then  the  incident  phasor  remains  unchanged  after  the  Lorentz  transformation. 
Multiplication  of  the  resultant  phasor  by  the  term  is  also  necessary. 

Lorentz  transformation  can  be  considered  also  as  rotation,  apart  from  multiplication  by 
,  but  then  axis  and  angle  of  such  rotation  depend  on  the  incident  polarization  point.  Axis  of 
that  rotation  is  perpendicular  to  the  great  circle  plane  through  the  I  and  T  points.  That  axis  is 
also  perpendicular  to  the  ‘imaginary  rotation’  axis  oriented  along  the  sphere  diameter  through 
the  I  point.  The  imaginary  rotation  axis  is  independent  of  the  T  point  location.  There  are  two 

rotation  matrices  corresponding  to  models  of  the  Jones  and  Sinclair  matrices:  Cg°^°  and 
* ,  accordingly. 

However,  there  is  a  possibility  to  use  conjugate  versions  of  those  matrices,  * 

and  ,  after  regular  inversion  and  Lorentz  transformations  in  other  Sinclair  and  Jones 
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matrices  (see  Appendix  C).  Necessity  of  using  those  other  scattering  and  propagation  matrices 
becomes  obvious  when  considering,  for  example,  the  regular  Jones  matrix  for  which  there  is  a 
need  to  find  its  special  polarization  points.  One  can  examine  the  other  Sinclair  matrix, 
,  for  the  same  obstacle  in  order  to  find  special  polarization  points  for  that  matrix 
using  well  known  methods  applicable  to  scattering  matrices.  Then  the  correspionding  special 
polarization  points  for  the  original  regular  Jones  matrix  can  be  immediately  found  by  applying 
to  that  ‘scattered’  phasor,  -  ,  the  orthogonality  transformation  according  to  the  rule: 

Ay,=c’^A^y^c^c^Ay, 

=  A"C’‘C’‘wf  =  =A''wf 

What  one  has  to  realize  is  that  special  polarization  points  of  the  A^  Sinclair  matrix  change 
their  meaning  after  orthogonalization  procedure.  Points  representing  incident  polarizations 
corresponding  to  maximum  emerging  power  exchange  their  location  with  those  producing 
minimum  output  power,  and  vice  versa.  Next,  co-pol-nulls  exchange  their  meaning  with 
cross-pol-nulls  what  means  that  in  the  case  of  propagation  the  cross-pol-nulls 
(eigenpolarizations)  always  do  exist,  as  co-pol-nulls  for  scattering,  while  the  polarizations 
which  become  orthogonal  may  sometimes  disappear. 

Interesting  is  the  case  of  the  scattering  matrix  model  in  the  so  called  characteristic 
basis  K.  Complex  elements  of  its  second  diagonal  are  of  opposite  sign  (their  sum  disappears) 
and  elements  of  the  first  diagonal  are  of  the  same  phase  (as  elements  of  symmetrical  matrix 
taking  diagonal  form  in  that  K  basis).  It  means  that  for  opix>site  direction  of  transmission  the 
model  rotates  by  180°  about  its  axis  (together  with  all  special  polarization  points).  The 
axis  of  rotation  after  inversion  is  always  situated  in  the  Vj,.  plane,  and  the  chord  joining 
co-pol-null  points  is  parallel  to  the  V,,  axis.  There  is  an  allowed  region  inside  the  Poincare 

sphere  model  in  which  the  inversion  I  point  can  be  situated.  Using  coordinates  of  the  I  point 
the  whole  Sinclair  matrix  in  the  AT  basis  (except  of  the  general  phase)  can  be  restored,  in  a  part 
of  that  allowed  region  two  solutions  exist  leading  to  two  different  Sinclair  matrices 
corresponding  to  the  same  I  point.  Depending  on  the  number  of  solution  and  on  location  of 
the  I  point  inside  the  allowed  region  different  mutual  locations  of  special  polarization  points 
result  Most  interesting  are  mutual  locations  of  those  points  for  the  inversion  I  point  situated 
on  the  boundary  surfaces  of  the  allowed  regions  for  the  two  solutions.  Special  geometrical 
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constructions  have  been  developed  enabling  one  to  find  special  polarization  points  for 
assumed  coordinates  of  the  I  point  located  on  boundary  surfaces  of  its  allowed  regions  or  in 
coordinate  planes  of  the  Poincare  sphere  corresponding  to  the  characteristic  K  basis  (see  [3]). 

The  nonsymmetrical  Sinclair  matrix  model  is  much  more  complicated  than  the  model 
of  symmetrical  scattering  matrix  considered  by  Kennaugh.  Kennaugh’s  model  has  dealt  with 
one  solution  only  for  the  scattering  matrix,  and  the  allowed  region  for  the  I  point  in  his  case 
was  only  the  negative  part  of  the  OQj^.  axis  including  the  center  of  the  sphere  (for  the  so 
called  ‘symmetrical’  targets)  and  the  point  on  its  surface  (for  the  ‘dipole  type’  targets). 

6.  CONCLUDING  REMARKS 

Presented  work  extends  the  Kennaugh’s  inversion  point  concept  to  the  case  of  Jones  and 
nonsymmetrical  Sinclair  matrices.  Polarization  and  phase  vectors  of  waves  interpreted  as 
tangential  polarization  phasors  have  been  used  to  construct  geometrical  models  of  those 
matrices  in  the  form  of  the  Poincare  spheres  of  the  same  diameter  and  equipped  with  common 
inversion  point  but  different  axes  and  angles  of  rotation.  Propagation  (optical)  and  scattering 
(radar)  transformations  using  Jones  and  Sinclair  matrices,  accordingly,  have  been  presented  in 
both  matrix  and  quatemionic  forms  exhibiting  full  agreement  of  the  two  approaches.  The 
choice  of  the  approach  can  be  left  to  researchers  who  will  apply  one  more  suitable  to  solve 
particular  problems. 
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APPENDIX  A.  COLUMN  PP  VECTOR  TRANSFORMATIONS 

Passive  TPP  sphere  rotation  transformation  -  change  of  the  PP  vector  (or  TP  phasor)  basis 


Denoting:  =[l,  ly\4,  and  = 


T 

'o' 

0 

,  U  — 

1 

,  one  obtains 


L  =u  and  I„  =  u  .  Therefore; 


where 


/.  /,k=k  ""k 
„■  «“k=l“"  «"k  k  «"k}' 
<]• 


4 


'X^*B 


w;  =ciu 


x''b 


a  -  b* 

rx  J“ 

- 1 

♦ 

_ 1 

Jx 


(A.1) 


(A.2) 


Active  rotation  transformation  -  change  of  the  PP  vector  (or  TP  phasor)  by  rotation  of  the 
TPP  sphere 


Indeed: 


Let 


if  C^.b=C^bC’ 

c^c‘«;  =c‘'<  =[»;■ 


^'P.A 


exp(-j^p  ) 
0 


0 

exp(+y^p  ) 


(A.3) 

(A.4) 


(A.5) 


then: 


~  ^  P,P  ~  ^  P^  P,A^  A 


a  -  b* 

A 

exp(-y^p  ) 

0 

a*  b* 

b  a* 

P 

0 

exp(+yVp  )_ 

-b  a 

aa  *  exp(-y^p  )  +  bb*  exp(+ j(/>f  ) 
a*  b  exp{-j(f>f  )-a*b  exp(+J(f>p  ) 

cos^f  -  Jqp  sin^f 

_  (Vp  -  .A>p)sin^p 


ab  *  exp(-j(f>p  )-ab*  exp(+j(f>p  ) 
-bb*  exp(-y^p  )  +  aa*  exp(+ j(/>f  ) 


A 


P 


(-Vp  -yu^)sin  (f>f 
cos^p  +  Jqp  sin 


(A.6) 


Also 


A  p  —  L.  p 


CA^k 

P.A^  A 


cosJ>p  Jq^  sin  ^p 


(-v^  -  Ju'*^)sm(f>^ 
cosipp  +  jq^  sin  (pp 


(A.7) 
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Therefore,  components  of  the  unit  axis  vector  for  the  real  rotation  are: 


A 

q 

A 

u 

= 

u 

= 

"2 

V 

P 

V 

K 

3- 

(A.8) 


Spatial  reversal  tran^ormation  (passive  or  active,  by  the  spatial  180°  rotation  about  the  y 
axis) 

Such  transformation,  passive  or  active,  is  of  the  same  form  and  can  be  presented  as 


^<=(CX)*  and/or 


When  applied  twice,  it  brings  the  vector  to  its  original  form. 


where 


with 


and 


Wg=Cg*C|Mg  and/or 


-w 

u 

B 

■-1 

o' 

U 

w  * 

_ ^0  _ 

—  ^x  ~~ 

X 

0 

1_ 

u  =  u®  =  {ab  *  +ba*Y^ 

=  sin  2Yx  cos2(^® 

w  =  w^=X-A^)^ 

=  cos  2^®  cos 2^®  cos2£:®  -  sin  26 ^  sin  26^ 

-  y(cos  2yx  26 x  sin  2ex  +  sin  26 x  cos  2ex  ) 


(A.9) 


(A.  10) 


(A.  11) 


(A.  12) 


(A.  13) 


Jones  and  Sinclair  matrices,  as  well  as  Lorentz  and  inversion  matrices,  both  represent 
active  transformations. 
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APPENDIX  B.  USEFUL  FORMULAE  FOR  RELATIONS  BETWEEN  SOME  SIMPLE 

AMPLITUDE  MATRICES 


For  the  following  matrices, 
C*  = 


0  -1 

1  0 


c°  = 


-iB 


-  w  u 


u  w 


;  detC°=-l,  = 


a  -h' 
b  o* 


satisfying  relations 


C*C’'  =C°*C°  *C^°^=diag{\,\} 

=-C^  C°  =C° 


;  detC""^  = 
(B.l) 


(B.2) 


easy  to  check  products  are  worth  noticing; 

^ROT  _ ^ROT  ^ 

^x^ROT  _ ^ROT 

^x ^ROT  ^  _ ^ROT ^x 


and 


C^C°*  =  C°C' 
C*C°  =C°  *C' 

c*c*  =  c°c* 


C’'C°C’'  =c*c°c*  =c°  * 

^  X  ^  O  ^  ^  X  ^  X  ^  ^  X  ^  o 
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Similarly,  for. 
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/ig  —  ^  /ig  , 


.INV  _  -X  .1 
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nB 
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nB 


in  any  B  basis  one  obtains: 


AT  +c^AT 


C*  =2 


1  0 
0  1 


and  A 


Tc’‘+c’‘aT*  = 


1  0 
0  1 


(B.6) 


Therefore,  from  definitions  of  the  Jones  and  Sinclair  matrices: 


As  =e^‘’ 


0  ^ROT  jLOR 
^B  ^nB 


and  Ag  =e-'^ 


_2_qROT  >|c  ^^OR 

j  '  B  nB 


(B.7) 


through  identities: 
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*C^  =  0^°^° C' A^^^  *C’ 


(B.8) 


and 


(v'')c“  ^cr‘-<Tc\  -^c'{a,  •e")=c’'crc‘'*=cr  'c'az  • 


(B.9) 

one  immediately  obtains. 

pL{(.4|e-'«-)+C’'(/i;  •e'")c>}=cr"  and  -)=iA,e^^^):'  aC'(a,  •<!'«)}=cr  * 

(BIO) 


APPENDIX  C.  DIFFERENT  JONES  AND  SINCLAIR  MATRICES  AND  THEIR 

MUTUAL  RELATIONS 


When  applying  relations  of  Appendix  B,  the  following  representations  of  different  Jones  and 
Sinclair  matrices,  useful  for  applications,  can  be  obtained; 


Original  and  modified  Jones  matrices: 


A° -pjf”  j. _ -c®  * /I  -r°*r*A^’  a°iF-2^u 


BO 

'b 


(C.l) 


A'^  - pji  2L _ 4^'^ -r^  A  -  r°  *  A”^  -  r°  * A°  a*iF  -2^ u 


sx 

‘b 


Original  and  modified  Sinclair  matrices: 

ClAl  =CtC^AT;  AX  =2^(ulr 

C^Al  ^ClAl  =C“CM,  =-^“;  (C.2) 

A’^uI=X^{uT)\  uT^-uT 

with  mutual  relations'  concerning  phase  factors, 

eJ^°  =-je^^-  (C.3) 
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rotation  matrices  (active  transformations), 


ROT 

B 


♦c*  =  jC’^ClC 
=  =  jC^C°  * 


(C.4) 


and  (active)  transformations  employing  the  inversion  I  point, 


.LOR  .m'.  aINV  a 


LOR 

nB 


(C.5) 


Rotation  quaternions: 


APPENDIX  D.  SOME  USEFUL  QUATERNIONIC  FORMULAE 


5 


B 

'X  ■ 

B  . 
'X 

C/-- 


--  cos  <f)]^ 


C, 


+  ^2^2  '*’^^3^3)2:  sin 
zo%<j)l  -(«,e,  +«2‘?2  sin^^ 

+("1^1  +"2^2  -”3^3)^^^  sin^^ 

_  ( yt  o  JL-  ft  o 


■■COS<f)x  +(«ie,  +«2^2 
:  COS(f>x  -  ('^i^i  +  "2^2 


^^3^3) A"  sin 


B 


(D.l) 


Lorentz  and  inversion  transformation  quaternions: 


=  1  -  y(Q^.  +  Ue2  +  We,)\ 
=  1  -  me,  +  Ue2  -  Ve3 ) 


I 

nB 


- “^3(1  “ y(Ql^i  +ulif2  +  v^V,)^} 
=h?"'‘^3 =ii-y(Q>.  +  u;e.  +v>,),)e, 
•=-<!,  (1  -  y(Q>,  +  u>.  -  v>, ), ) 
'i'"'  •  =  HT  *  =  (1  -  Xq;«,  +  u;.!.  -  V>,).  )e, 


(D.2) 


(D3) 


Stokes  quaternions  (of  incident  waves): 
for  use  with  Jones  matrices, 

^B  =^r  + A(ie,  +ue^  +ve3)^}<t>  V2  =  ^ 


1  +  q  u  -  yv 
u  +  /V  1  -  q 


(D.4) 


for  use  with  Sinclair  matrices, 

PI’  ^Pl*  =  ji{^+J((\e,  +\xe^-ve^)l}<:x>42nl  *ul  =~ 


■  1 

1  +  q  u  +  /V 
u  -  /V  1  -  q 


(D.5) 


Stokes  quaternions  of  emerging  waves  are  of  the  type  for  use  with  Jones  matrices 
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APPENDIX  E.  MUELLER  AND  KENNAUGH  MATRICES  -  A  COMPARISON 
In  the  elliptical  basis  B  of  collinear  co-phase  phasors  the  two  matrices  take  forms, 
Kennaugh  matrix  (see  Appendix  A  for  C®  matrix): 


K„  = 


a, 

biB 

^3B 

'l 

1 

0 

0  ' 

C,B 

azB 

^4B 

=  U(Ab<S>Ab*)U  =  D^^K^; 

u=-j= 

0 

0 

I 

-y 

^3B 

C4B 

^3B 

^2B 

V2 

0 

0 

I 

j 

C6B 

^2B 

^4B  ^ 

I 

-I 

0 

0  _ 

D^=U*(C^®C|*)U 


(E.I) 

(E.2) 


For  Ag  = 


A, 


and  M *,  ^L2,3,4,  the  resulting  formulae  for  elements  of  the 


Kennaugh  matrix  are: 


a,  +M,  +M,  +M,),  b„  -M,  +M,  -M,), 

a,,  -M,  -M,  +M,),  c„  +M3  -M,  -M,), 

asB  =Re(A^A^*+A^A^*)g 
=Re(^4^3 

bjB  =Im(>l2>I,  *+A^A^*)b 
CjB  =lm(A^A^  *-A^A^*)b 

bjB  =Re(>l2^3  *+A,A^*)B  c^B  =Re(A,A^  *+AjA,*)b 

^4B  =  Re(^2^3  *  -^4/1, *)g  C4B  =  Re(^2^4  * -AjAj *)b 

b,B=  ImiA^Aj  *  +A^A, *)b  c,b  =  ImiA^A,,  *  +AjA, *)b 

b,B  ^ImiA^A,  *-A^Aj*)b  c^b  =lm(A2A^  *-A,A,*)b 


(E3) 


(E.4) 


(E.5) 


Mueller  matrix  (see  Appendix  A  for  Cg  matrix): 

b ,  D  b 


K°b  = 


1 

o 
IB 
O 

3B 

A 


IB 

o 

2B 

o 

4B 

o 

6B 


'3B 

b:^ 

^  O 


3B 

O 

2B 


D°b=U(C°<8>C°*)U  = 


'SB 

Kb 


'2B 


=  u*(A!<s>a^*)u=diKb 


0  0 

l-2u^  -2uRew 
-2uRew  u^-Rew^ 


0 

-  2u  Im  w 
-  Im 


-2u  Im  w 


Im  w 


u 


Re 


(E.6) 


(E.7) 


Z.H.  Czyz  ONR-Report-7 (Final),  November  23,  2002,  N00014-02-1-0222,  02PR04940-00  40 


Elements  of  the  first  raw  of  Kennaugh  and  Mueller  matrices  are  identical. 


For  Al  = 


^2  ^3  ^  A=l, 2,3,4,  the  resulting  formulae  for 

A°  A° 


elements  of  the  Mueller  matrix  are: 


a,  +M°  +m:  -M°  +m: -M°), 

a2B  =i(A^2°  -A^3°  -K  +K)b  +^3° -M: -M°), 

a°,  =Ke(A:A°  *+A°Ar)s 
a:,=Re(A:Ar-A°Ar)s 

b°,=Im(A:A°*+A°A°% 


c°,=ImK^3°*-^°A°*)B 


b3,=Re(^3X*+^4X*)B 
b:=Re(^3X*-^4X*)B 
b3,=Im(^3X*+^4X*)B 
b°  =Im(A°A°  *-A°A°*)g 


cl=Re(A°A:*+A°Ar), 

c:,^Re(A°A:*-A°Ar)s 

c°,=Im(A°A:*+A°A°*), 

c:,=Im(A°A:*-A°A,°*), 


(E.8) 


(E.9) 


(E.IO) 
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APPENDIX  F.  DETERMINATION  OF  SPECIAL  POLARIZATION  POINTS  FOR 

SCATTERING  MATRIX 


Efficient  method  determines  polarization  ratios  of  special  incident  polarization  points 
(Compare  [11]); 


= 


NT„Ta 


R. 


expi 


};  A  =  R,"  +  RjR, 


The  corresponding  Stokes  parameters  are: 


'q' 

TpTj 

co%2y 

u 

= 

sin  2y  cos  2(5 

V 

B 

sin  2y  sin  25 

(F.l) 


(F,2) 


Parameters  determining  special  polarization  points: 

M  and  N=Mx  (incident  polarizations  for  maximum  and  minimum  emergent  power) 

R,  =i(A,A,  *-A,A,  *+A,A,*-A,A,*),  =b„ 

=-(^,4  *+4^2*)fl  =-R3*  =  -b3fl  +J^5B, 


Oi  and  O2  (co-pol-nulls). 


•^1  ~  -^3fl  -^4fl 

~^A^g,  =  — 2^2fl  5 


(F.4) 


El  and  Ez  (cross-pol-nulls  or  eigenpolarizations);  they  exist  if  A  >  0 , 

R,={A,A,*-A,A,*+A,A,*-A,A,*)g 

R,=2{A,A,*-vA,A,*)g  (F.5) 

R,=2{A,Ar-^A,A*)g 


K  and  L=Kx  (co-pol-max/saddle  points  or  characteristic  polarizations  determining  the  Qj^. 
axis), 

r,=(a,a,*-aa*)b 

R,^-A,g(A,*+A,*)g-A,g*(As+A,)g^R,*. 


In  the  case  of  polarizations  M  and  N  simpler  formula  can  be  more  convenient  directly 
expressing  their  normalized  Stokes  parameters: 


(F.7) 
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